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Abstract 

It is shown that the problem of designing a two-reflector system 
transforming a plane wave front with given intensity into an output 
plane front with prescribed output intensity can be formulated and 
solved as the Monge-Kantorovich mass transfer problem^. 

1 Introduction 

Consider a two-reflector system of configuration shown schematically on Fig. 

1. Let ( X — [xi, X2, Xn), z) be the Cartesian coordinates in M'^"''"'^, n > 

2, with z being the horizontal axis and Xi,X2, ■■■,Xn the coordinates in the 
hyperplane a : z = 0. Let Bi denote a beam of parallel light rays propagating 
in the positive direction and let Q denote the wavefront which is the cross 
section of Bi by hyperplane a. Assume that f2 is a bounded domain on a. An 
individual ray of the front is labeled by a point x G 0. The light intensity of 
the beam Bi is denoted by I{x),x G Q, where / is a non- negative integrable 
function. 
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Figure 1: Sketch for Problem I 



The incoming beam Bi is intercepted by the first reflector Ri, defined 
as a graph of a function z{x),x G Cl. The rays in Bi are reflected off Ri 
forming a beam of rays B2. The beam B2 is intercepted by reflector R2 
which transforms it into the output beam ^3. The beam ^3 also consists 
of parallel light rays propagating in the same direction as Bi. The output 
wavefront at a distance d > from the hyperplane a is denoted by 7^; the 
projection of Ta on the hyperplane a we denote by T. The second reflector R2 
is also assumed to be a graph of a function w{p),p G T. The quantity 

where Pd is the map of fl on and J is the Jacobian, is the expansion ratio 
and it measures the expansion of a tube of rays due to the two reflections 
[5] . It is assumed that both Ri and R2 are perfect reflectors and no energy is 
lost in the transformation process. Consequently, the corresponding relation 
between the input intensity I on fl and output intensity L on is given by 

L{Pd(x))\J(Pd(x)\^I(x). (1) 

The "two-reflector" problem that needs to be solved by designers of opti- 
cal systems consists in determining the reflectors Ri and R2 so that all of the 
properties of the two-reflector system above hold for prescribed in advance 
domains fl,T and positive integrable functions I{x),x e fl, and L{p),p e T; 
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see Malyak [6] and other references there. It is usually assumed in applica- 
tions that Q and arc bounded and convex. 

Two fundamental principles of geometrical optics are used to describe 
the transformation of the beam Bi into beam B^: the classical reflection 
law leading to the ray tracing equations defining the map Pd and the energy 
conservation law for the energy flux along infinitesimally small tubes of rays; 
see [6] where the problem is formulated for rotationally symmetric data and 
a class of rotationally symmetric solutions is found. 

The problem of recovering reflectors Ri and R2 without assuming rota- 
tional symmetry was formulated rigorously by Oliker and Prussner in [7] , and 
it was shown that it can be considered as a problem of determining a special 
map of n — T with a potential satisfying an equation of Monge-Ampere 
type relating the input and output intensities. Existence and uniqueness of 
weak solutions were established by Oliker at that time but only the numer- 
ical results implementing a constructive scheme for proving existence were 
presented in [7] for several test cases. Detailed proofs were given in [8]. 

In this paper we show that this problem can also be studied in the frame- 
work of the Monge-Kantorovich mass transfer problem studied by Brenier [2] , 
Caffarelli [3], Gangbo and McCann [4], and other authors. In our notation, 
the Monge-Kantorovich mass transfer problem is to transfer the intensity / 
on Q into the intensity L on T via a map P : Q — > T for which the total 
transportation cost f^C{x, P{x)) Idx is minimized. Here C{x,p) is a given 
strictly convex cost function. 

The proof of existence and uniqueness of solutions to the Monge- 
Kantorovich problem is obtained by solving a minimization problem for the 
functional 



considered on pairs of continuous functions on and 00 onT that satisfy 



Under various conditions it is shown in [2], [3], [4] that this functional is 
minimized by some pair (Co)<^o) (referred to as Kantorovich potentials), and 
that P(x) = X + VCo solves the Monge-Kantorovich problem. 

Applying these ideas, we show that the geometric optics problem at hand 
can be formulated as a Monge-Kantorovich mass transfer problem with a 





C{x) - uj{p) > -C{x,p), X eQ, peT. 



(3) 
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quadratic cost function; see section 6. The Kantorovich potentials corre- 
spond to the pair of reflectors that solve the problem. The condition (3) has 
a geometric meaning; namely, it filters out reflectors that allow only optical 
paths longer than a certain prescribed one. The functional (2) to be mini- 
mized is the mean horizontal distance between points of the two reflectors, 
with the average weighted by the two intensities. 

We prove that there are always two different reflectors system satisfying 
the stated requirements. The corresponding ways in which one intensity is 
transfered into the other one are exactly the most and the least energy efficient 
in the sense of the Monge-Kantorovich cost. This result is thus ultimately a 
variant of Fermat's principle. 

The fact that the solution to the above geometrical optics problem can 
be derived from a variational principle gives rise to a numerical treatment 
of the problem different from the one suggested in [8]. In particular, when 
the problem of minimizing (2) under constraints (3) is discretized we have 
the linear programming problem with quadratic constraints. We ran some 
numerical experiments with this approach and intend to return to this point 
in a separate publication. 

This paper is organized as follows. In section 2, we recall some results 
from [8] concerning the ray tracing map, assuming smoothness of the reflec- 
tors, and formulate the main "two-reflector" problem. In section 3, we give 
a geometric characterization of reflectors as envelopes of certain families of 
paraboloids. Such characterization is of independent interest. In section 4 
we use this geometric characterization to define weak solutions of type A and 
type B of the two-reflector problem. To prove existence and uniqueness of 
solutions for each type we utilize the ideas of the Monge-Kantorovich the- 
ory and introduce the functional (2) on a certain class of "quasi-reflector" 
systems. This is done in section 5. In the same section it is shown that the 
problem of finding weak solutions of type A is equivalent to finding minimiz- 
ers of (2). Weak solutions of type B correspond to maximizers of (2). On the 
other hand, existence of minimizers (maximizers) to this functional is not 
difficult and has been established before in [2], [3], [4]. This implies existence 
of solutions. The uniqueness in respective class is established in section 6 
by proving that the ray tracing map P associated with a weak solution min- 
imizes or maximizes the quadratic Monge-Kantorovich cost for which the 
functional (2) is the dual. The main theorem on existence and uniqueness of 
weak solutions to the two-reflector problem is stated and proved in section 
7. 
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Finally, we note that same methods lead to variational formulations and 
solutions of other geometrical optics problems involving systems with single 
and multiple reflectors. 



2 Statement of the problem 

We begin by reviewing briefly the analytic formulation of the problem for 
smooth reflectors; see [8] for more details. 

Let Ri be given by the position vector ri{x) — {x,z{x)), x & Cl, with 
z e C^{fl). The unit normal u on Ri is given by 



Consider a ray labeled by x E fl and propagating in the positive direction k 
of the z— axis. According to the reflection law the direction of the ray y{x) 
reflected off Ri is given by 

y — k — 2{k,u)u — k — 2- ^ ' ^ 



l + \Vz\ 



where (, ) is the inner product in 7?"+^. Denote by t{x) the distance from 
reflector Ri to reflector R2 along the ray reflected in the direction y{x) and let 
s{x) be the distance from R2 to the wavefront along the corresponding ray 
reflected off i?2- Assume for now that t e C^{Q) and R2 is a hypersurface. 
The total optical path length (OPL) corresponding to the ray associated 
with the point x E Cl is l{x) = z{x) + t{x) + s{x). A calculation shows that 
l{x) = const = I on fl. Since 

R2 : r2{x) = ri{x) +t{x)y{x), x e Cl, (4) 

the image of x on the reflected wavefront is given by 

Pd{x) = ri{x) + t{x)y{x) + s{x)k, x e n. (5) 

The equation (5) is the ray tracing equation for this two-reflector system. 

Introduce the map P{x) — Pd{x) — dk : Q ^ T. A calculation [8] shows 
that 

p ^ P{x) ^ X + l3Vz{x), xefl, (6) 
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where (3 = I — d is the "reduced" optical path length. 

To simphfy the notation we will write L{P(x)) instead of L{Pd{x)){= 
L{P{x) +dk)). For the input intensity /(a;), x eVL, and the output intensity 
L[P[x)) on Td we have in accordance with the differential form of the energy 
conservation law (1), 

L{P{x))\J{P{x))\^I{x), xen, (7) 

where we also take into account that J{Pd) — J{P)- It follows from (7) that 
D,, T, I and L must satisfy the necessary condition 

/ L{p)dp^ [ I{x)dx. (8) 
Jt Jn 

It follows from (6) that J{P) = det [Id + PHess{z)], where Id is the 
identity matrix and Hess is the Hessian. Hence, by (7), 

L{x + (3Vz)\det[Id + f3Hess{z)]\^I in O. (9) 

Thus, the problem of determining the reflectors Ri and R2 with properties 
described in the introduction requires solving the following 

Problem I. Given bounded domains fl and T on the hyperplane a and two 
nonnegative, integrable functions I on D, and L on T satisfying (8), it is 
required to find a function z e C^{Cl) such that the map 

Pa^x + (5Vz : ^ f (10) 

is a diffeomorphism satisfying equation (9). 

It is shown in [7] that once such a function z is found, the function w 
describing the second reflector is determined by z and (5 as 

w{P{x)) = d - s{x) = z{x) + -dV^I' - 1). (11) 
Following [8] we introduce the function 
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Then by (6) and (11) 
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V - {x,VV) + hwV\^ 



P = VV, (13) 

(14) 



where V — {x, W) is the negative of the usual Legendre transform of V. 
Thus, y is a potential for the map P : — > T. If P is a diffeomorphism then 
the inverse of the transformation {x,z{x)) — > {p,w{p)), where p = P(a;), is 
given by 

x{p) = p-\p) =p- (3Vpw{p), (15) 
z{p)^w{p)-^\VMp)\' + ^), P^T. (16) 
In terms of the potential V the equation (9) becomes 

L(VV)\detHess(V)\ ^ I in Q, (17) 

which is an equation of Monge- Ampere type. 

In order to clarify the relations between the parameters I, d and P, note 
first that it is the reduced optical path length /3 that is intrinsic to the 
problem. The choice of cross sections of the fronts (that is, the selection of 
a particular value for d) is extraneous. In fact, it is easy to see that if {z,w) 
are two reflectors as above and we change d to d' then I' — d' = l — d = (3 
and {z, w) are not affected by such change. 

Finally, we note that the two-reflector system described above has the 
following two symmetries. For A e put z'{x) — Xz{x) and P' — jP. Then 

P'{x) = P{x), 
w\p) = Aw(p) + ^(A- i), 

where P'{x) — x + (3'Vz\x). In other words, the system is invariant un- 
der some combination of flattening (stretching) the flrst and translating and 

flattening (stretching) the second reflector. 

Note also that a horizontal translation of both reflectors, that is, adding 
the same constant to z and w, does not change /3 and the map P. 
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3 Geometric characterization of reflectors 



We examine first more closely the relationship between the functions V 
and w for smooth reflectors. For the rest of the paper we assume that and 
T are hounded domains on the hyperplane a. We continue to assume that 
the map P : Q — > T is a diffeomorphism. Let x e p e T and 

Q{x,p) = {x,p) + (5w{p) - —. 

If p = P{x) then by (12) - (14) we have 

V{x)^Q{x,P{x)). 

Denote by Sy the graph of V over fi. Let G and po — P{xq). The 
tangent hyperplane to Sy at (xq, ^^(^^o)) is given by the equation 

Z = {x,po) - {xo,Po) + V{xo), 

where {x, Z) denotes an arbitrary point on that hyperplane. Taking into 
account (14), we obtain 

(x,po) - {xo,Po) + V{xo) = {x,po) + pw{po) - y = Q{x,po). (18) 

Since V{xo) = Q{xo,po), we conclude that Z = Q{x,po) is the tangent 
hyperplane to Sy at {xo,V{xo)). Consequently, if V is convex then Z — 
Q{x, po) is a supporting hyperplane to Sy from below and if V is concave then 
Z = Q{x,Pq) is supporting to Sy from above (relative to positive direction of 
the 2;-axis). Because T is bounded there are no vertical tangent hyperplanes 
to the graph of V and we have 

V{x) > Q{x,po) for all x e Q iiV is convex, 
V{x) < Q{x,po) for all X E iiV is concave. 

Since for every p & T the hyperplane Q{x',p) is supporting to Sy at some 
(x', V{x')), we get 

V{x) > Q{x,p) for all X E Cl, p E T ii V is convex, (19) 
V{x) < Q{x,p) for all X E fl, p E T if is concave (20) 
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and in both cases we have equahties Hp — P{x). 
Let 

U{p) = — - I3w{p) - y, R{x,p) = {x,p) - I3z{x) - —. 

It follows from (19) and (20) that 

U{p) > R{x,p) for all a: G p e T if is convex, 
U{p) < R{x,p) for all a; e 0, p e T if is concave 



(21) 
(22) 



and in both cases equalities are achieved if p = P{x)- Also, for any fixed 
Xq E fl the hyperplane R{xo,p) is supporting to the graph Su of U at (po = 

P(xo),f/(po)). 

Using the usual characterization of convex functions [10] we obtain from 
(19), (21) and (20), (22) 

V{x) — sup Q{x,p), U{p) = sup R{x,p) when V is convex, 

peT xen 

V{x) — mf_Q{x,p), U{p) — inf R{x,p) when V is concave. 
For convex V this implies 

z{x) — sup 
p&f 



w(p) = inf 



(5"^ — \x — 
\x — pp — 13"^ 



+ w{p) 
+ z{x) 



X eQ. 
pef. 



(23) 
(24) 



Similarly, when V is concave we have 



z{x) — inf 
per 

w{p) = sup 



- \x -p\ 
\x — p\ 



+ w{p) 



2/3 



+ z{x) 



X eil. 
peT. 



(25) 
(26) 



The characterizations (23), (24) and (25), (26) have a simple geometric 
meaning. To describe it, consider first the case when V is convex. Recall 
that the total optical path length I — z{x) + t{x) + d — w{P{x)) — const (see 
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Fig. 1). Also, t'^{x) = \x-P{x)\'^ + \z{x)-w{P{x)\^ and p = l-d. It follows 
from (23) that {x, z{x)) is a point on the graph of the paraboloid 

f3^ — \x — pP 

kp,w{x) = — h X ea, (27) 

with the focus at {p = P{x),w{P{x))) and focal parameter /?. 

Similarly, it follows from (24) that a point (p = P{x),w{P{x))) on the 
second reflector lies on a paraboloid 

\x — pP — 

hxAp) = ^ pea, (28) 

with the focus at {x, z{x)) and focal parameter /3. 

Let be the convex body bounded by the graph of paraboloid 

kp,w{x) and Hx,z{x) the convex body bounded by the graph of paraboloid 
hx,z{p)- Then (23) and (24) mean that the graphs of z{x) and of w{p) 
are given by 

Sz^d(^\j K,,^^p^ , (29) 

= y j . (30) 

When the potential V is concave wc have similar characterizations of ^S";;; 
and with |J in (29), (30) replaced by f]. 

Remark 3.1. It follows from (23), (24) that when V{x) is convex then for 
any x e fl the path taken by the light ray through the reflector system is 
the shortest among all possible paths (not necessarily satisfying the reflection 
law) that go from {x, 0) to {x, z{x)), then to some (p, w(p)) and then to {p, d). 
Of course, the shortest path satisfles the reflection law and p — P{x). For 
concave V the corresponding light path is the longest as it follows from (25), 
(26). Thus the characterizations (23), (24) and (25), (26) are variants of the 
Fermat principle. 



4 Weak solutions of Problem I 

We use the geometric characterizations of reflectors in section 3 to define 
weak solutions to Problem I. Let fl and T be two bounded domains on the 
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hyperplane a and (5 a fixed positive number. 

Definition 4.1. A pair [z, w) E C(0) x C(T) is called a two-reflector of type 
A if 

z{x) = sup/cp,^(p)(x), X efl, (31) 
per 

w{p) = inf p e f , (32) 

where A;p „,(p)(a;) and hx,z{x){p) are defined by (27) and (28). Similarly, a pair 
{z, w) e C{Vt) X C(r) is called a two-reflector of type B if 

2;(x) = inf fcp_^(p)(x), X e fi, (33) 

w{p) = sup/i^,^(^)(p), pef. (34) 

To avoid repetitions, we consider below only two-refiectors of type A. 
The changes that need to be made to deal with two-refiectors of type B are 
straight forward and arc omitted. 

It will be convenient to construct the following extensions z* of the func- 
tion z and w* of w to the entire a. For a pair {z,w) as in definition 4.1 
let 

y(x) = y + /3^(x)-^ and g(x,p) = (x,p)+«p)-^. (35) 

It follows from (31) that 

V{x) — sVi^Q{x,p), x e 

That is, V is convex and continuous over Vl. Furthermore, since T is bounded, 
the graph Sy has no vertical supporting hyperplanes. For any fixed p & T 
define the half-space 

Q+{p) = {{x, Z)eaxR' I Z> Q{x,p)}. 

Then 

Sv*^dlf]Q+{p) 
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is a graph of a convex function V* defined for all x E a. Note that V*{x) — 
V{x) when x & fl. We now define an extension of z by putting 



z%x) 



Similarly, for any fixed x E fl we let 

R+(x) = {(p, Z) eaxR^ I Z> R{x,p)}, 



where 



X 



R{x,p) = {x,p) - pz{x) - —, pea. 



Then the function 

U*{p) = sup R{x,p), pea, 
is defined. It is also convex. It follows from (32) that for p e T 

U*{p)^^-Pw{p)-^ {^U{p)). 
The corresponding extension of w we define as 



w*{p) 



»2 



pea. 



Lemma 4.2. The function V* is uniformly Lipschitz on a with Lipschitz 
constant maxy Also, U*{p), p e a, is uniformly Lipschitz on a with 
Lipschitz constant maxfj In addition, z e Lip{Cl) and w e Lip{T) with 
the Lipschitz constant < ^"P(^.p)egx^ '^"^l . 

Proof. By our convention the normal vector to a plane Q{x,p) (when p is 
fixed) is given by (— p, 1). It follows from definition of V* that Sy has 
no supporting hyperplanes with normal {—p, 1) such that p ^ T. Since T 
is bounded, this implies the first statement of the lemma. The statements 
regarding U* are established by similar arguments. From these properties of 
V* it follows that the function z* is continuous on a and Lipschitz on any 
compact subset of a. Similar properties hold also for w*. 

Now we estimate the Lipschitz constant for z on Q. Let {z, w) be a two- 
reflector of type A. Let x,x' e and let z{x') > z{x). (If the opposite 
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inequality holds we relabel x and x' .) Fix some small e > Q. It follows from 
(31), (32) that there exists a p' e T such that z{x') < kp/^yj(j,'){x') + e. Then 

\z{x') - z{x)\ < kp,^^(^p,){x') - z{x) + e < kp,^^(^p,){x') - kp, ^^(^^p,){x) + e 

< sup \V kpi ^yj{^i){x)\\x' — x\-\- e = — sup \s — p'\\x' — x\ + e 
x&ci P sen 

< — sup \s — p\\x — x\ + e. 

Letting e — > 0, we obtain the statement regarding Lipschitz constant for 
z. The same statement regarding w and two-reflectors of type B are proved 
similarly. □ 

Next, we define the analogue of the ray tracing map P for a two- 
reflector. For that we need to recall the notion of the normal ma,p [1], p. 114. 
Let M : G — be an arbitrary convex function defined on domain G G a 
and Su its graph. For xq & G let 

Z - u{xo) = (p, X - Xo) 

be a hyperplane with normal (— 1) supporting to Su at {xq,u{xq)). The 
normal map : G — > a at is defined as 

^u(xo) = IjM, 

where the union is taken over all hyperplanes supporting to Su at {xq, u{xq)). 
Definition 4.3. Let {z, w) be a two refiector of type A. For x E a we put 

P{x) = z/y*(a;). 

For reflectors of type B the ray tracing map is deflned similarly with the 
use of function —U*. In general, P may be multivalued. 

Lemma 4.4. Let {z, w) be a two-reflector of type A and z* and w* their 
respective extensions as above. Then P{x) G T for all x E a. In addition, 
for any p E T the set {x E fl\P{x) = p} ^ 0. Furthermore, for any x E fl 

P{x) = {all pef \ w{p) = /i.,.(.)(p)}. (36) 
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Proof. Let x & a and Q{x,p) a supporting hyperplane to V* at {x,V*{x)). 
Then the normal p is in z/y. (x). On the other hand, by definition of V*, Sy* 
has only supporting hyperplanes with normals in T. Hence, P{x) C T. 

Let p E T and Q{x,p) a supporting hyperplane to S'y*. We need to show 
that there is an x e Q such that p e P{x). By (35) and (32) we have for any 
X eO. 

V{x) - Q{x,p) = P (/i.,.(.)(p) - w{p)) > 0. 

By (32) there exists a,n x e Cl such that V{x) — Q{x,p) — 0. This implies 
the remaining two statements of the lemma. □ 

Remark 4.5. It follows from definition that P is multivalued at points x 

where Sy* has more than one supporting hyperplane. At such x the function 
z* is not differentiable. Let (xq, 2*(a;o)) be one such point. Then 

P{xq) — {p E a \ Q{x,p) is supporting to Sy* at {xo,V*{xo))}. 

In other words, a light ray labeled hj Xq E Q that hits a point where the first 
refiector has a singular point will split into a cone of light rays. These rays 
will generate a subset on the paraboloid hxQ,z{xo){p) whose projection on a 
is P{xo). This is consistent with the physical interpretation of diffraction at 
singularities of this type [5] . 

Remcirk 4.6. Since V* is convex, by Rademacher's theorem, the Lcbesgue 
measure of the set of singular points on Sy* is zero. Thus, P{x) is single- 
valued almost everywhere in a. Furthermore, the functions z (z*) is a dif- 
ference of two convex functions, and therefore, it is differentiable almost 
everywhere in (a). The same is true for w and w*. It follows then from 
the definitions of P and V that for almost ell x E Q 

P{x) = VV = x + pVz{x), (37) 

A similar property holds also for the function w. 

Lemma 4.7. If {z, w) is a two-reflector of type A then for all x E Cl, p E T 

z{x)-w{p)>j^i(3'-\x-p\'). (38) 

In addition, for almost all x E Q there exists a unique p E T such that 
p = P{x) and (38) in this case is an equality. 
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Proof. The lemma follows from (31), (32), Remark 4.6 and Lemma 4.4. □ 
Define the inverse of P for p G T as 

P~^{p) = {x e a I p e P{x)}. 

Theorem 4.8. Let B be the a-algebra of Borel subsets of T. Let {z, w) be 

a two-refllector of type A. For any set t E B the set P~^{t) is measurable 
relative to the standard Lebesgue measure on a. In addition, for any non- 
negative locally integrable function I on a the function 

C{t) = / I{x)dx 
Jp-Ht) 

is a non-negative completely additive measure on B. 

Proof. The proof of this theorem is completely analogous to the proofs of 
Theorems 9 and 16 in [9]. □ 

Lemma 4.9. Let and T be two bounded domains on a and I a non- 
negative integrable function on Q expended to entire a by setting I{x) = 
for X E a\VL. Let {z,w) be a two-reflector of type A or B. Then for any 
continuous function h onT we have the following change of variable formula 

[ h{p)C{dp) = I h{P{x))I{x)dx. (39) 
JT Jn 

Proof. In the integral on the right h{P{x)) is discontinuous only where P is 
not single valued, that is, on the set of measure zero. Thus, the integral on 
the right is well defined. 

We may assume that J^I{x)dx > 0; otherwise, the statement is triv- 
ial. Fix some small £ > and a positive integer N. Partition the interval 
[min/i(p), max/i(p)] into sub-intervals (Si, (Sjv of length < e/ J^I{x)dx and 
let hi e Si. Put Tj = {p e T \h{p) e Si}. Then for sufficiently large N 

I j^h{p)C{dp)-Y,hiC{Ti) \<e. 
1, iV, i ^ j, meas{P-\Ti) f]P-\Tj)) = (see Remark 4.6). 
I / h{P{x))I{x)dx -y^hi / I{x)dx\<s. 



For any i,j = 
Hence, 
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This, together with the previous inequahty, imply 

I / h{P{x))I{x)dx-y2hi^ri)\<2e. 
Jn 

□ 

Definition 4.10. A two-reflector {z,w) of type A ( B) is called a weak 
solution of type A (B) of the two-reflector problem I if the map P : Q — > T 
is onto and for any Borel set r C T 

C{t) = j^L{p)dp. 

Using Lemma 4.9 and this definition we obtain 

Lemma 4.11. Let {z,w) be a weak solution of type A (B) of the two-reflector 
problem I. Then for any continuous function h on T 

I h{p)L{p)dp= I h{P{x))I{x)dx. (40) 

5 A variational problem and weak solutions 
of the two-reflector problem 

As before, we consider here only the case of two-reflectors of type A. We 
comment on the case of two-reflectors of type B in section 7. Let, as before, 
/ and d be the given parameters of the system, and P = I — d. Let ( G C{fl) 
and (jj e C{T). With any such pair (C,<^) we associate a "quasi- reflector" 
system in which the light path is defined as follows. Let {x,p) & Cl x T and 
let {x, C{x)) be the point where the horizontal ray emanating from {x, 0) hits 
the graph of (. Let P2 be the point where the horizontal ray begins in order 
to terminate in (p, d). The "light" path is the polygon [x, 0) P1P2 {p, d); see 
Fig. 2. The length of this path is 

l{(,uj,x,p) = ({x) + V {({x) - uj{p)y + \x-p\'^ + d- uj{p). 

The class of admissible pairs is defined as 

Adm(n,T) ^{(Cu) eC(n) xC{f) I l(C,u;,x,p)>l y(x,p) eClxf}. 

(41) 
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Figure 2: Definition of / {(,u!,x,p). 

By construction, a pair {C,uj) G C(fi) x C{T) lies in Adm{Q,T) if and 
only if for all x G p G T 

C(x) > A;p,^(p)(x), < K,i:{x){p), (42) 

or, equivalently, if and only if 

C{x)-u;{p)>^{f3'-\x-pn (43) 

for all x e Q, p e T. It follows from (31), (32) that a two-refiector of type A 
is an admissible pair. 

The following functional is central to our investigation. 

Let / and L be two nonnegative integrable functions on Q and T, respec- 
tively, satisfying the energy conservation law (8). For (Ceo) e Adm{0,,T) 
put 

J^{C,uj)^ / C{x)Idx- / uj{p)Ldp. 
Jn Jt 

Clearly, is linear and bounded on Ci^Cl) x C(T) with respect to the 
norm max{| |^| |oo, | Ico"! |oo}- Geometrically, J-'{Cy '-^) is proportional to the mean 
horizontal distance between the points of the two graphs, the average being 
weighted by the intensities. 

We consider now the following 
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Problem II. Minimize T on Adm{Q,T). 

Proposition 5.1. [2, 3, 4j There exist {z,w) e Adm{^,T) such that 

Tiz.w) — inf TiC.oj). 

(C,w)GAdm(n,T) 

We may further assume that the pair (2;, w) satisfies the conditions (31), (32) 
for a type A reflector system. 

Proof. It follows from (42) that one can restrict the search for minimizers to 
two-reflectors (C,'^) of type A. Also, because of (8 is invariant under 

translations C, C, + p, uj ^ uj -\- p ior constant p e R. Hence, it is sufficient 
to consider only two-reflectors (C, oS) for which C,{xq) — for some xq G Ct. 

By Lemma 4.2, C, and uo are uniformly Lipschitz with the Lipschitz con- 
stant K = supj.gQpgy \x — p\/ p. It follows that for all x e O 

|C(x)| = |C(x)-CMI<i^diama 

For all p G T we have 

1 1 

uj{p) < /ixo,C(xo)(p) = ^(ko <max — (|xo-?|^-/3^). 

Finally, since /ix,c(a;)(p) = (k - P? - /3^)/2/3 + C{x) > -(3/2 - Kdiamfi for 
all a; e 0, we also get the lower bound 

x&l Z 

Therefore, the maps C and cu are also uniformly bounded. By the Arzela- 
Ascoli theorem, and because J-' is continuous, the infimum of J- is achieved 
at some two-reflector pair {z,w). □ 

Wc now show that weak solutions of Problem I and sohitions of the Prob- 
lem II are the same. This establishes existence of weak solutions to Problem I. 

Theorem 5.2. Let z e C(Q), w e C{T) be a two-reflector of type A. Then 
the following statements are equivalent. 

(i) {z,w) minimizes T in Adm{Q,,T). 

(a) {z,w) is a weak solution of type A of the two-refiector Problem I. 
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Proof, (ii)^(i). Let {C,uj) e Adm{n,T). Then by (43) and Lemma 4.7 for 
almost all X e n 

Cix) - ujiP{x)) >^(^p^-\x- P{x) |2) = z{x) - w{P{x)). 

Integrating this inequality we get 

/ C{x)Idx - / uj{P{x))Idx > / z{x)Idx - / w{P{x))Idx. 
Jn Jn Jn Jn 

Using Lemma 4.11, we get 



J-'{(,uj) = / ({x)Idx — / uj{p)Ldp> / z{x)Idx — I w{p)Ldx = J-'{z,w). 
Jn Jt Jn Jt 

Since (C, t^) was arbitrary, we are done. 

(i)^(ii). Let P denote the ray tracing map for the pair (-2, tf). It must 
be shown that Jp_T^^^^ I{x)dx = J^L{p)dp for all Borel sets t CT. We will 
prove that this is the Euler-Lagrange equation for the functional J-. 

It is sufficient to establish this for the case when r is an open ball with 
center po e T and radius r > contained in T. For i — 1, 2, • • • define for 
p E a 

{1, if\p-po\<r -\ 

i(r-|p-po|), iir - \ <\p-po\ <r 
0, if |p -pol > r. 

Then Xi is continuous on a, < Xi ^ Ij and the sequence {Xi{p) • P G 
converges on a pointwise to the characteristic function of r Xt{p) ■ 
Fix some i and for e e (—1, 1) put 

w,(p) ^ w(p) + £ ■ Xi(p) 

Ze{x) = sup kp,y,^(p){x) = sup{— {Po-Ix-pI"^) +w,{p)}. 
peT peT 

By construction, the pair {z^, w^) satisfies the condition (43), with ^ replaced 
by z^ and cu repaced by ty^. We show now that z^ belongs to Lip(fi). Let 
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x,x' & Q and let z{x') > z{x). (If the opposite inequality holds we relabel x 
and x'.) Let p' e T be such that z^{x') — kp/^ws(p'){x')- Then 

^ kp' ,we(p'){^ ) ~ kpi ^yj^(^pi-j(x) 

< sup \^kp',w^{p'){s)\\x' — x\ = -5 sup \s — p'\\x' — x\ 
sen P sen 

< — sup |s — — x\. 
P sen,peT 

Hence, is continuous and (z^jW^) G Adm{fl,T). 

Now let a: G n. For each e let be a point in T such that z^{x) — 
kpe,we(pe)i^) ■ (Thls cholcc, of course, may not be unique.) Then 

Ze{x) — Z{x) = kp^,We{pe)i^) " ^(^) — kp^,ws{ps)i^) ~ kp^^u)(ps){^) 

Similarly, if p G P{x), then 

^e{^) ~ ^{^) — ^e{^) ~ kp^u>{p){^) ^ kp^u}e{p){^) ~ kp^y}(jp){^x) 

= £Xi(p)- 

Therefore, 

-\e\ < exiip) < Ze{x) - z{x) < exiipe) < kl (44) 

for all X G fi. In particular, z^ converges uniformly to 2; on Q as £ —> 0. 

Now consider those a; G for which the ray tracing map P is single- 
valued. This is the case for almost all x G fi. We claim that p^ = P^, where 
Ps denotes the ray tracing map for {zg,Ws), converge to p as £ — > 0. 

Suppose this is not true. Then there is a sequence {Pej}, j = 1,2,... 
with £ j ^ as J — > 00 and a constant 77 > such that 

|p — Psj \ > 1] for all j. 

Let z'{x) = maXp/g2^ |p,_p|>^ /cp'^^(p')(a;). Note that the maximum in the 
definition of z' is attained. Since p is the unique point in T such that 
z{x) — kp^^(j)){x), it follows that z' < z{x). Therefore, for all j 

z{x) - z,.{x) = z{x) - A;p,^,,^,^.(p,^.)(a;) = z{x) - A:p,.,«,(p,.)(a;) + w{p,.) - w,.(p^.) 

= Z{x) - A;p,^,,„,(p,.)(x) - SjXiiPej) 
> z{x) — z! — |£j|. 
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This contradicts to the fact that Zs.{x) converges to z{x) as j — > oo. Hence, 
^ p if £ ^ 0. 
It follows from (44) that 

z^{x) - z{x) 



- Xi{p) < \Xi{Pe) - Xi{p)\- 

Letting £ — > and using the continuity of Xi, we conclude that for almost all 

X efi 

Ze{x) = Xi{p) = Xi{P{x)). 

6=0 



de 



Then 



d_ 

de 



Ze{x)I{x)dx = I Xi{P{'^))I{'^)dx. 

e=0 J ^ Jii 



Since T has a minimum at {z, w), we obtain 
d 







de 



J^{ze,w^)^ / Xi{P{x))H^)d^ ~ / Xi{p)L{p)dp. 

£=0 Jt 



Now let i — >• oo in this equality. This is possible as Xiip) ~^ Xr{p) pointwise 
on a and therefore Xi(-P(a;)) — > x,-(P(x)) pointwise almost everywhere on fi. 
Then, noting that for almost all x e Q = Xp-i(t)(^)' we obtain 

/ L{p)dp= / Xr{p)L{p)dp^ / Xr{P{x))I{x)dx 
Jt Jt Jq 

= / Xp-'^(T)i^)H^)dx = / I{x)dx. 
Jn Jp-'^{r) 

This completes the proof of the theorem. □ 



6 Connection between the two-reflector and 
Monge-Kantorovich problems 

In our notation, the Monge-Kantorovich mass transfer problem [2] can be 
formulated as follows. Consider the class of maps P : Q T which are 
measure-preserving, that is, they satisfy the substitution rule 

[ h{P{x))Idx= I h{p)Ldp 
Jn Jt 

for all continuous functions h on T. Each such map is called a plan. 
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Problem III. Minimize the quadratic transportation cost 



P^l [ \x - P{x)\'^Idx. (45) 

among all plans P. 

Note that for weak solutions to the two-reflector Problem I the ray tracing 
map P is a plan by Lemma 4.11. In fact, we have the following 

Theorem 6.1. Let {z,w) be a weak solution of type A of the two-reflector 

Problem I. Let P be the corresponding ray tracing map. Then P minimizes 
the quadratic transportation cost (45) among all plans P : Q ^ T, and any 
other minimizer is equal to P almost everywhere on supp(/)\{a; e fl : I{x) = 
0}. 

Proof. Let P : fl ^ T he any plan. Then by (38), 

z{x) - w{P{x)) > ^{/3' -\x- P{x)n (46) 

for almost ell x E Q and equahty holds if and only if P{x) — P{x). Integrat- 
ing against Idx and applying Lemma 4.11, we get 



J_ I [fi^-\x-P{x)\^]Idx< [ [zix) - w{P{x))]Idx 

z{x)Idx — / w{p)Ldp = / [z{x) — w{P{x))]Idx 
i Jt Jn 

[f3^ -\x- P{x)\^]Idx. 



j 

Jq 



n 

_ 1 

This shows that P is a minimizer of the transportation cost. 

To show uniqueness, note that if equality holds in the integral inequal- 
ity, then equality must hold in (46) for almost all x e supp (/) \ {/ = 0}. 
Therefore, P = P a.e. on supp(/) \ {/ = 0}. □ 

Remark 6.2. The functional !F is the dual of the quadratic cost functional 
(45) [2]. 
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7 Existence and uniqueness of weak solutions 
to the two-reflector problem 

Theorem 7.1. There exist weak solutions of type A and of type B to the 
two-refiector Problem I. Furthermore, if [z, w) and {z', w') are two solutions 
of the same type with ray tracing maps P and P' , respectively, then 

P{x) = P'{x) 

for almost all x e supp(/) \ {x e Q : /(x) = 0}. 

Proof. By Proposition 5.1 and Theorem 5.2 we know that Problem I has 
a solution. The only property that remains to be checked is that for that 
solution, P : Q — > T is onto. But this follows from Lemma 4.4. 

In order to prove uniqueness, let [z, w) and {z' , w') be two solutions of 
type A, with ray tracing maps P and P' , respectively. Then by Proposition 
6.1, both are minimizers of the quadratic Monge-Kantorovich cost functional, 
so that P — P' a.e. on supp(/) \ {/ = 0}. 

So far, existence and uniqueness of weak solutions to the two-refiector 
Problem I has been shown for weak solutions of type A. However, it is clear 
that a similar result holds for weak solutions of type B as well. Namely, 
we can define A(im+(0, T) as the space of all pairs of refiectors such that 
l{^,ui,x,p) is less than /, and then that T admits a maximizing pair on 
Adm+{Vt,T), which is a weak solution of the two- reflector problem I. This 
shows, in particular, that for such a solution the ray tracing map maximizes 
the quadratic transportation cost among all plans. □ 

Corollary 7.2. Suppose that in addition to the assumptions in Theorem 7.1 
the function I > m Q. Then there is a constant p e R such that 

z'{x) = z{x) + p 
w\p) = w{p) + p 

for all X & Q and all p E T . In other words, weak solutions of type A are 
unique on VL and T up to a translation of the reflector system, and the same 
result holds for type B solutions. 

Proof. We show this for weak solutions of type A. By the theorem, P{x) = 
P'{x) for all X & fl. Note that by (37), Vz{x) — Vz'{x) for almost all x E fl. 
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It follows that there is a constant p such that z'{x) — z{x) + p for all x e Jl. 
Now, by definition of w'{p), 



w'{p) = inf /ia:,z'(x)(p) = inf K,z{x){p) + p = w{p) + p 



for aU peT. 



□ 
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